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Ice formation is one of the most common and important processes on Earth and almost always occurs at 
the surface of a material. A basic understanding of how the physicochemical properties of a material’s 
surface affects its ability to form ice has remained elusive. Here we use molecular dynamics simulations to 
directly probe heterogeneous ice nucleation at an hexagonal surface of a nanoparticle of varying hydrophilicity. 
Surprisingly, we find that structurally identical surfaces can both inhibit and promote ice formation and 
analogous to a chemical catalyst, it is found that an optimal interaction between the surface and the water 
exists for promoting ice nucleation. We use our microscopic understanding of the mechanism to design a 
modified surface in silico with enhanced ice nucleating ability. 


Upon cooling, liquid water crystallizes into solid ice. 
Due to the presence of a free energy barrier separating 
the liquid and crystalline states, however, it is possible 
for liquid water to remain in a metastable ‘supercooled’ 
state to temperatures far below the equilibrium melting 
temperature. Heterogeneous ice nucleation, that is, ice 
nucleation in the presence of impurity particles such as 
mineral dust, soot or certain types of bacteria, gener¬ 
ally increases the rate of ice nucleation and is the domi¬ 
nant process by which ice forms in nature. 1 Recent work 
has argued that any ice formation at temperatures above 
—20° C must necessarily occur heterogeneously.^ Empir¬ 
ically, a large variance in the propensity of different ma¬ 
terials to nucleate ice is observed, and due to the impor¬ 
tance of ice formation in e.g. the climate sciences, much 
effort has been expended in identifying and cataloging 
the effectiveness of different materials to nucleate iceP 
This has motivated many simulation studies of heteroge¬ 
neous ice nucleation in the pr esen ce of different surfaces, 
including graphite kao'.iniu ’ ^and silver iodide.® De¬ 
spite the vast amount of research into heterogeneous ice 
nucleation, major gaps in our knowledge still exist, espe¬ 
cially with regard to our understanding of the underlying 
chemical physics; this is reflected in our inability to accu¬ 
rately predict a material’s ice nucleating efficiency and to 
answer seemingly simple questions such as how does hy¬ 
drophilicity affect the ice nucleation rate ? Not only is an 
understanding of the chemical physics of heterogeneous 
ice nucleation needed to pr edic t the ice nucleating effi¬ 
ciency of existing materialsp®^ but it is also paramount 
for the rational design of new materials to either pro¬ 
mote or inhibit ice nucleation. Controlling ice formation 
is desirable in a variety of fields, for example, in the cry- 
opreservation of cells, tissues and organspP the food and 
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transport indu stries and even as a potential means for 
climate control. 12 13 

In contrast to fields such as chemical catalysis 1 ^ and 
materials design}- 15 there is currently no comprehensive 
set of design principles in terms of molecular ‘descrip¬ 
tors’ for making new substances to control ice forma¬ 
tion. Put more simply, we do not know which are the 
relevant microscopic properties of a material that deter¬ 
mine its macroscopic ice nucleating efficiency. Often, the 
so-called ‘requirements’ for a good ice nucleating agent 
(INA) have been discussed, such as the requirement for 
a good crystallographic match to ice and the ability of 
water to chemically bond to the surface of the particle 
(i.e. hydrophilicity) P3 Although properties such as a 
good crystallographic match are impo rtant in heteroge¬ 
neous nucleation of some systems} 1 A 1 8 I such criteria have 
neither served as a full set of guidelines to identify good 
INAs, 1 nor have they aided the systematic improvement 
of ice nucleation inhibitors or promoters. Experimentally 
there is disagreement regarding the role of hydrophilicity. 
For example, Alizadeh et aZ™ have found ice nucleation 
to be slower on superhydrophobic surfaces, which they 
attribute not only to a lower contact area between the 
water and the surface, but also to a larger free energy 
barrier to nucleation. In contrast, Li et a/ ! 20 l 2l l found ice 
nucleation to be enhanced at hydrophobic modified sili¬ 
con wafers relative to their unmodified hydrophilic coun¬ 
terparts, which was attributed to a faster dynamics of 
water at the hydrophobic interface. Recently, Lupi et al. 
investigated the role of hydrophiclity o f gr aphitic surfaces 
using molecular dynamics simulation: 3 by varying the 
hydrophilicity in different ways (by uniformly changing 
the interaction of water with the surface or by introduc¬ 
ing hydrophilic species at the surface), they found that 
the ice nucleating efficacy of the surface could either in¬ 
crease of decrease. Also, it is found on kaolinite (a known 
hydrophilic INA) and platinunJ® that the most stable 
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water ove rlaver can inhibit the growth of subsequent wa¬ 
ter layersP^El Furthermore, in the case of requiring a 
good crystallographic match, evidence for ice-like struc¬ 
tures at surfaces is in general lackingpH and there are also 
instances where materials with a good crystallographic 
match to ice are poor INA^SET ( we no t e that recent sim¬ 
ulation studies 7 ® have found the unreconstructed basal 
face of silver iodide to act as a template for ice). We are 
therefore either faced with the prospect of relying on ex¬ 
periments to determine the efficacy of IN As on a case-by- 
case basis, or we can try and rationalize their behavior 
by elucidating the underlying molecular processes that 
control heterogeneous ice nucleation. 

Here, in the first of a series of two articles, we present 
results from molecular dynamics (MD) simulations where 
we directly probe heterogeneous ice nucleation in the 
presence of a face centered cubic (FCC) nanoparticle 
(NP). The NP exposes its hexagonal (111) surface as its 
principal facet and can therefore act as a template for the 
hexagonal basal face of ice. With this NP completely im¬ 
mersed in water, shown in Fig. [lja), we perform a series 
of studies in which we systematically explore the depen¬ 
dence of the nucleation rate on the hydrophilicity of the 
NP. By comparing these results to reference simulations 
of homogeneous nucleation we find a very interesting de¬ 
pendence of the nucleation rate on NP hydrophilicity; 
the NP can both promote and inhibit ice nucleation and 
exhibits a maximum nucleation rate at intermediate in¬ 
teraction strengths with the water. By examining the 
molecular level details of the nucleation processes at dif¬ 
ferent hydropliilicities of the NP, we find that the struc¬ 
ture in the immediate vicinity of the interface couples 
strongly with the nucleation rate. We then use this un¬ 
derstanding of the underlying chemical physics to design 
an improved INA. This first article emphasizes how we 
can use our microscopic understanding to control ice nu¬ 
cleation. In the second article, we discuss certain aspects 
of the mechanism in greater detail, as well as contextual¬ 
izing this work with respect to previous studies on surface 
hydrophilicity and ice nucleation. 

We have used the single site mW potential to model 
the interactions between water molecules,®^ which al¬ 
lows us to investigate length- and time-scales inacces¬ 
sible to ice nucleation simulations that employ more tra¬ 
ditional empirical potentials.® The NP was modeled as 
an FCC crystal with a lattice constant of 0.392 nm, con¬ 
sisting of 380 atoms. Previous work has suggested that 
such a lattice may aid in structuring water into ice-like 
arrangements!^! The FCC NP was hemispherical and ex¬ 
posed its (111) face as its primary facet (approximately 
2.5nm in diameter). For the interaction between the NP 
atoms and the water molecules, the Lennard-Jones po¬ 
tential U{r) = 4e [(<r/r) 12 — (cr/r) 6 ] was used, where r is 
the distance between an NP atom and a water molecule. 
The hydrophilicity of the NP was controlled by varying 
e (a constant value of a = 0.234 nm was used through¬ 
out). As mentioned earlier, Lupi et aZpBI controlled the 
hydrophilicity of graphitic surfaces not only in this man¬ 
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FIG. 1. (color online) (a) Snapshot of a typical ice nucleation 
event on the NP. Ice-like molecules are colored blue and the 
NP is colored silver. The NP is totally immersed in water 
(liquid-like molecules are shown by gray dots), (b) Varia¬ 
tion of the nucleation rate with the strength of the water- 
NP interaction J5 a d s - As E^ds increases, so too does the hy¬ 
drophilicity. The solid blue line indicates homogeneous nu¬ 
cleation (the dotted lines are an error estimate): data above 
and below this line indicates promotion and inhibition of ice 
nucleation by the NP, respectively. It can be seen that at 
weak and strong water-NP interaction strengths the NP in¬ 
hibits nucleation, while at intermediate interaction strengths 
(A a ds/AHvap ~ 0.15-0.6), the NP strongly promotes nucle¬ 
ation. 


ner, but also by introducing hydrophilic species at the 
surface, and found opposite trends: we discuss the pos¬ 
sible causes of this apparent discrepancy in more detail 
in our second paperP 2 ! Interactions were truncated after 
0.753 nm. This setup yielded contact layers at a height 
between 0.2-0.25 nm above the (111) surface of the NP, 
which is in reasonable agreement with values obtained 
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FIG. 2. (color online) Sensitivity of the water structure at the surface of the NP to the water-NP interaction strength, (a) Ice 
nucleation at the NP with E^ds/AH vap « 0.3. The water molecules in contact with the NP (colored red) form an hexagonal layer 
commensurate with the surface that resembles the basal face of ice. The water molecules directly above this contact layer (colored 
blue) also form a similar hexagonal structure, (b) The structure of water at the NP with E a d B /AH vap sa 0.08. The water-surface 
interaction is too weak to stabilize an ice-like structure, (c) The structure of water at the NP with _E a d s / AH vap « 1.2. The 
water surface interaction is too strong and the water molecules cannot rearrange into an ice-like configuration. In both (b) and 
(c), nucleation occurs away from the surface in a homogeneous manner. 


from density functional theory calculations of water at 
metal surfaces G233I1 We must emphasize, however, that 
we are using simplified model surfaces in order to under¬ 
stand possible general trends that may underlie heteroge¬ 
neous ice nucleation and that one must exercise caution 
in trying to make one-to-one correspondences with actual 
surfaces. 

All simulations were performed using the LAMMPS sim¬ 
ulation packagd® with 2944 mW molecules in a peri¬ 
odic supercell. Previous simulation studies have sug¬ 
gested that the critical ice nucleus varies from ~10 water 
molecules at 180 K® to ~85-265 at 220 fj—IAI giving us 
confidence that our simulations should not be subject to 
serious finite size effects. Furthermore, simulations using 
a slab geometry (approximate dimensions of 66 A 2 ) with 
4000 mW molecules confirm that the conclusions drawn 
from this work are not affected by changes to the box 
size and shape. For each value of the water-NP interac¬ 
tion energy, 16 MD simulations were performed at 205 K 
and 1 bar. Under these conditions, bulk liquid mW water 
is still metastable (as opposed to unstable but under¬ 
goes homogeneous nucleation on a timescale accessible 
to computer simulation such that statistically meaning¬ 
ful rates can be obtained. To detect ‘ice-like’ molecules, 
we have used the CHILL algorithm of Moore et a/.® (As 
the CHILL algorithm was designed for bulk homogeneous 
nucleation, it does not necessarily capture the full behav¬ 
ior in regions of broken symmetry i.e. interfaces. Never¬ 
theless, it is useful as a qualitative visual aid.) By mon¬ 
itoring the potential energy, we are able to determine 
the induction time to nucleation for each simulation and 
thus the probability Pn q (t) that a given system remains 
liquid after a time t from the start of the simulation. 


We are able to determine the ice nucleation rate R by 
fitting -Piiq(t) = exp [—(Jtt) 7 ], where 7 > 0 is also a fit¬ 
ting parameter. Further details of the fitting procedure 
and simulation setup are provided in the Supplemental 
Material.® In order to gauge the effectiveness of the NP 
as an INA, we have also studied bulk homogeneous nu¬ 
cleation, using identical settings. 

Explicit simulations of heterogeneous ice nucleation 
have only recently started to emerge in the literature 
(see e.g. Refs. IM30I1 and[3Hl) and to enable a system¬ 
atic study, we draw conclusions from over 200 success¬ 
ful nucleation trajectories in total. Fig. [ljb) shows the 
dependence of the nucleation rate on the water surface 
interaction, the main finding of this study. Specifically, 
we have plotted log 10 (I?/I?hom) vs E a d s /Aff vap , where 
-Rhom is the bulk homogeneous rate and A H vap is the en¬ 
thalpy of vaporization of bulk mW water (10.65 kcal/mol 
at 298 K)A rich variety in the ice nucleating behavior 
is seen: the NP is seen to both promote and, surpris¬ 
ingly, inhibit ice formation. We expect this inhibition 
effect to be concentration dependent; as the NP con¬ 
centration becomes more dilute, we expect the rates to 
tend to that of homogeneous nucleation. At low val¬ 
ues of E a( j s , the heterogeneous nucleation rate is approx¬ 
imately two times lower than i?hom- Thus when the par¬ 
ticle is very hydrophobic, it tends to inhibit nucleation. 
As the water-surface interaction strength increases so too 
does the nucleation rate until it reaches a maximum at 
Uads/Af/vap ss 0.4 that is nearly 25 times faster than 
bulk homogeneous nucleation. Beyond the maximum, 
the rate steadily decreases until £ a d s /AH vap « 1.0. Fur¬ 
ther beyond this, the rate remains roughly constant and 
slightly below i?hom- 
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We now try to understand this intriguing dependence 
of the nucleation rate on the hydrophilicity of the NP. 
To this end, we have examined in detail the mechanisms 
by which nucleation occurred on the NP for the vari¬ 
ous interaction strengths. As the (111) surface of an 
FCC crystal exhibits hexagonal symmetry, one possible 
mechanism for heterogeneous ice nucleation is a tem¬ 
plate effect whereby the molecules in the contact layer 
form an hexagonal structure commensurate with the sur¬ 
face. Fig. §a) confirms this, where we show a typi¬ 
cal ice nucleation event in the presence of the NP with 
Sads/AiJvap ss 0.3 (close to the maximum rate). Here we 
can clearly see that the water molecules in contact with 
the (111) surface of the NP do indeed form an hexagonal 
structure commensurate with the surface that resembles 
the basal face of ice. We can also see that the water 
molecules directly above this contact layer also form a 
similar hexagonal structure. The surface is therefore act¬ 
ing to promote ice nucleation by providing an arrange¬ 
ment of adsorption sites that resemble the structure of 
ice, thereby stabilizing structural fluctuations towards 
ice-like arrangements in the liquid. 

Now that we have established that the NP acts to pro¬ 
mote ice nucleation by acting as a template for ice, the 
dependence of the rate on E a d s can easily be understood 
as a competition between water-water and water-surface 
interactions. In Fig. [2j)b) we show the structure of wa¬ 
ter in contact with the NP for E a d s /AH vap ps 0.08 (the 
weakest E a d s investigated, which inhibits ice nucleation). 
Clearly, such a weak water-surface interaction is unable 
to stabilize ice-like configurations and in fact, ice nucle¬ 
ation is seen to occur away from the surface in a homoge¬ 
neous manner. Willard and Chandler have found that the 
structure of the interface between water and a hydropho¬ 
bic substrate is akin to the liquid-vapor interface® a re¬ 
cent simulation study from Haji-Akbari et al. 41 has found 
that ice nucleation is disfavored at the liquid-vapor inter¬ 
face. This appears to be consistent with our observations 
and with those of Lupi et al. in Ref. [3] Fig. 0c), on the 
other hand, shows the structure of water in contact with 
the NP for E a d s /AH vap « 1.2. While this NP also in¬ 
hibits ice nucleation, it does so for the opposite reason: 
the water-surface interaction is too strong, meaning that 
the water molecules cannot rearrange to form an ice-like 
layer at the surface. It is also clear that the coverage 
is higher than when ice forms at the (111) surface, as 
shown in Fig. 0a) (we also show this quantitatively in 
the second paper in this serietf®). For this strongly inter¬ 
acting scenario, we also see that ice forms away from the 
surface in a homogeneous manner. This is also consistent 
with the observations of Reinhardt and Doy^H on ice-like 
surfaces. The observed coupling between the molecular 
mechanism and the ice nucleation rate as we change the 
surface hydrophilicity is actually rather simple; we now 
demonstrate how we can exploit such simplicity to design 
a surface with improved ice nucleating efficiency. 

When an ice-like hexagonal overlayer forms at the 
(111) surface of the NP, such as in Fig. 0 a), it does so 



FIG. 3. (color online) Surface modification to promote 
ice nucleation for strong water-NP interaction strengths 
(A a ds/AHv a p > 1.0). By introducing small adsorbates (col¬ 
ored yellow) at the (111) surface, the template effect can be 
recovered (c.f. Fig. 0c)). The nucleation rate is increased 
by approximately a factor 50 compared to the bare NP for 
Aads/AHvap 1.2, as indicated in Fig. 0b)- 


with sub-monolayer coverage i.e. not all of the available 
adsorption sites are occupied by water molecules. We 
refer to these unoccupied adsorption sites on the (111) 
terrace as “excess” sites. For E a( ± s /AH vap > 0.6, these 
excess sites are occupied for long times and for nucleation 
to occur, an area of decreased coverage at the surface 
must occur such that an hexagonal motif can form. This 
motif can then act as a template for the hexagonal basal 
face of ice (a movie showing this for E a d s /AH vap ss 0.9 
is provided 3 -). When E a d s > AH vap it becomes favor¬ 
able for a water molecule to occupy a site on the surface, 
including the excess sites, rather than a position in the 
bulk liquid. This prevents the water molecules in the 
contact layer from forming the hexagonal arrangements 
required for ice nucleation at this NP. By this rationale, if 
the density of available adsorption sites was lower, then 
the template effect (and the enhanced nucleation rate) 
may be preserved at higher values of E a d s . To this end, 
we have modified the (111) surface of the NP by ad¬ 
sorbing small molecules, at the excess sites, which only 
have a weak interaction with water,® and recomputed 
the nucleation rate with E a d s /AH vap « 1.2. As seen in 
Fig- 0 the template effect is indeed recovered. It is also 
seen in Fig. 0b) that this modified surface enhances ice 
nucleation by a factor of 50 compared to the unmodified 
surface. This is a clear demonstration of how the molecu¬ 
lar insight into heterogeneous ice nucleation can be used 
to rationally design surfaces of different ice nucleating 
ability. Experimentally, this could be realized through 
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adsorption of small molecules to the surface (e.g. carbon 
monoxide) or through surface alloying. In fact, alloying 
a platinum (111) surface with tin is observed to promote 
the formation of an hexa gonal ice-like bilayer under ultra- 
high vacuum conditions j 44 l 45 l in a fashion analogous to 
our modified surface (the tin atoms at the surface act in 
part to reduce the density of adsorption sites). We also 
note that in our simulations, the nucleation rate can also 
be decreased by adsorbing small molecules such that the 
NP can no longer act as a template. 

More generally, the sensitivity of the nucleation rate 
on surface hydrophilicity could be tested by e.g. using 
nanoparticles of gold or silica functionalized with organic 
molecules of varying lrydrophobicity. In addition to us¬ 
ing well established methods such as the droplet freezing 
techniques,® it may also be possible to exploit recent ad¬ 
vances in femtosecond X-ray scattering techniques that 
have allowed real-time monitoring of homogeneous ice nu¬ 
cleation in micron sized water droplets.® Not only could 
such an experimental protocol be used to compare rates 
of ice nucleation in the presence of immersed NPs, but 
information regarding the impact of such NPs on the mi¬ 
croscopic structure of the liquid should also be available. 

In summary, we have used computer simulations to sys¬ 
tematically compare heterogeneous ice nucleation rates 
in the presence of a simple model nanoparticle of vary¬ 
ing hydrophilicity. This complements a number of recent 
simulation studies on specific systems.® 5 ® We have 
seen that the nanoparticle can promote ice nucleation 
by acting as a template for the hexagonal ice lattice, but 
that the ice nucleating efficiency is lost if adsorption is too 
strong, due to a high coverage of water molecules destroy¬ 
ing the template effect. Modification of the surface such 
that the coverage of water molecules is reduced recov¬ 
ers this template effect and enhanced nucleation can be 
achieved for strongly adsorbing surfaces, clearly demon¬ 
strating how molecular level understanding of heteroge¬ 
neous ice nucleation can be used to manipulate the rate of 
ice formation. The use of molecular descriptors to predict 
useful macroscopic properties of materials has been suc¬ 
cessfully used in other fields, such as chemical catalysis.® 
Designing new catalysts for reactions such as methana- 
tion (CO + 3H 2 —» CH 4 + H 2 0) has relied upon the 
establishment of a Sabatier principle based on a compu¬ 
tationally tractable quantity (in this case the dissocia¬ 
tion energy of CO at the surface).® We have seen that 
for the surface investigated in this study, the adsorption 
energy of a single water molecule can be used to describe 
the heterogeneous ice nucleation rate. Although a com¬ 
prehensive set of rules still requires further experimental 
and theoretical investigation, the results presented here 
suggest that if the surface acts as a template for ice, 
then one must tune either the density of adsorption sites, 
or the propensity of water to adsorb to the surface. In 
our second article,® we show that the variation of the 
ice nucleation rate upon surface hydrophilicity is depen¬ 
dent upon the surface topography, demonstrating that 
the combined effect of different surface properties needs 


to be considered when trying to understand what makes a 
good INA. Other properties such as the crystallographic 
match to ice and the role of surface defects are also likely 
to be important, as will more complex interactions such 
as electrostatics and explicit hydrogen bonding. The re¬ 
sults presented in this letter serve as a platform upon 
which future studies can be conducted. 
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